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This report is an attempt to explain both s- and p- 
wave nonleptonic hyperon decays by means of the QCD 
enhanced effective weak Hamiltonian supplemented by the 
SU(3) Skyrme model used to estimate nonperturbative matrix 
elements. The model has only one free parameter, namely, 
the Skyrme charge e, which is fixed through the experimental 
values of the octet-decuplet mass splitting A and the axial 
coupling constant gA- Such a dynamical approach produces 
nonleptonic hyperon decay amplitudes that agree with exper- 
imental data reasonably well. 

PACS number(s): 12.39. Dc, 12.39.Fe, 13.30.Eg 



In the Skyrme model, baryons emerge as soliton con- 
figurations of pseudoscalar mesons 0- ||. Extension 
of the model to the strange sector, in order to account 
for a large strange quark mass, requires that appropri- 
ate chiral symmetry breaking terms are included. The 
resulting effective Lagrangian can be treated by starting 
from a flavor symmetric formulation in which the kaon 
fields arise from rigid rotations of the classical pion field 
HHI?! . The associated collective coordinates are canoni- 
cally quantized to generate states that possess quantum 
numbers of the physical strange baryons f|||,^| ■ It turns 
out that the resulting collective Hamiltonian can be di- 
agonalized exactly even in the presence of the flavor sym- 
metry breaking (SB) jl. This approach leads to a good 
description of hyperon masses, charge radii, magnetic 
moments, etc. ||. It should be noted that in the first 
phenomenological applications of the Skyrme model one 
attempted to fit absolute baryon masses, which required 
a ridiculously small pion decay constant . Nowadays 
it is understood that there exist 1/N C corrections to the 
total baryon masses that are not fully under control and 
therefore only mass splittings can be reliably reproduced. 
In this approach f w is kept at its experimental value. 
Hence the results for nonleptonic hyperon decays (NHD) 
H need to be updated accordingly. 

Both s- and p-wave NHD amplitudes were quite suc- 
cessfully predicted bv using quark models with QCD en- 
hancement factors 



11 1. Note that there are not only 



current-algebra and ground-state exchange pole-diagram 
terms, but there exist other important contributions to 
both s and p waves. The so-called factorizable contribu- 
tions and/or kaon poles were estimated in Pole- 



diagram contributions to p waves from the (l/2 + ) - 
Rooper type of resonances and to s-waves through the 
(1/2 - ) -resonance exchange were calculated in [ p"2"[ . 

This report is an attempt to test whether the effective 
weak Hamiltonian and the extended SU(3) Skyrme model 
are able to predict both s- and p-wave NHD amplitudes. 
The minimal number of couplings Skyrme model is used 
to estimate only the nonperturbative matrix elements of 
the 4-quark operators |8[|. All remaining quantities en- 
tering the expressions for the decay amplitudes like mass 
differences, coupling constants etc. are taken from ex- 
periment. This approach uses only one free parameter, 
i.e., the Skyrme charge e. In order to avoid the unnec- 
essary numerical burden, throughout this report we use 
the arctan ansatz for the Skyrme profile function [T^ ] . 

The starting point of our analysis of NHD in the frame- 
work of the Standard Model || is the effective weak 
Hamiltonian in the form of the current ® current inter- 
action, enhanced by QCD. It is obtained by integrating 
out heavy-quark and VF-boson fields. This Hamiltonian 
contains the 4-quark operators Oi and the well-known 
Wilson coefficients [^,0. For the most recent values, see 
Ref. [|nj. For the purpose of this paper, we use the Wil- 
son coefficients from Ref. [[To): c\ = —1.90 — 0.61C, C2 = 
0.14 + 0.020C, c 3 = c 4 /5, c 4 = 0.49 + 0.005C, with 
C = V* d V ts /V* d V us . Without QCD short-distance cor- 
rections, the Wilson coefficients would be ci = — 1, C2 = 
1/5, C3 = 2/15, and C4 = 2/3. In this paper we simply 
consider both possibilities and compare the results. 

Note that there exists a different approach of Ref. 
|^4] , ^5| in which meson-baryon couplings are directly ob- 
tained from the chiral Lagrangian. There, the effective 
phenomenological constants extracted from experiment 
take into account all QCD effects hidden in the structure 
of the effective Hamiltonian (including the enhancement 
factors embodied in the values of the a constants) . This 
approach gives comparable results for the s- waves but 
fails for the p- waves |0| . 

The techniques used to describe NHD (1/2+ -> 1/2+ + 
_ reactions) are known as a modified current-algebra 
(CA) approach. The general form is 

<n( q )B'(p')\H e Jf\B(p)>=u(p')[A(q) 



<B'(p')\H w \B(p)> 



V(q) 



j 5 B(q)]u(p) 
-«%)• (1) 



Here the first term is the CA contribution, the second 
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is the modified pole term, and the third is a term that 
vanishes in the soft-meson limit. The V(q) term contains 
the contribution from the surface term, the soft-meson 
Born-term contraction, and the baryon pole term, which 
are combined in a well-known way J^,[w| . It represents a 
continuation of the CA result from the soft-meson limit. 
Further continuation is contained in the factorizable term 
S(q), which is proportional to the meson four- momenta. 

The parity-violating amplitudes A receive contribu- 
tions A c from CA commutator terms, factorizable terms 
S(q), and pole terms from the (l/2~) - resonance ex- 
change. The main contributions to the B amplitudes 
come from the baryon pole terms V(q), including both 
the ground state and the radially excited states. 

The current-algebra A c and baryon-pole B v ampli- 
tudes are well known from the literature. They contain 
weak matrix elements denned as <ibb' =< B'\H^ C \B >, 
which have the following general structure: 



a BB , = V2G F V: d V us < B'\ Cl Of c \B > 



(2) 



The factorizable term S(q) is calculated by inserting 
vacuum states. It is therefore a factorized product of 
two current matrix elements, where the first two-quark 
current is sandwiched between baryon states, while the 
second two-quark current is responsible for pion emission. 

The CA and the baryon-pole terms contain the 4- 
quark operator matrix elements, which are nonpertur- 
bative quantities. This is exactly the point at which 
the Skyrme model can be used. Each of the operators 
Oi from (^) contains four types of operators, namely, 
duus , dsuu , dsdd, dsss, and takes the form of the prod- 
uct of two Noether SU(3) currents, which can be found 
in Refs. f|,[l(|. In our calculations we use four operators 
Oi . The first of them is 

0\ = -<Zz,7^(Ai - iA 2 )<?L qLl^i^+ihhL, (3) 

where the SU(3) properties are expressed explicitly in 
terms of the Gell-Mann A-matrices. The connection with 
the effective Hamiltonian operators Oi is obvious. 

In order to estimate the matrix elements enterin g (|2| ), 
we take the SU(3) extended Skyrme Lagrangian p|p.6|: 
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(2) 
SK 



SB 



C 



wz, 



(4) 



d 4 xTr (d^Ud^) 



etc. 



(1 2) 

where Cg' K , £wz, and Csb denote the a- model, 
Skyrme, symmetry breaking (SB), and Wess-Zumino 
(WZ) terms, respectively. For U(x) £ SU{2), the SB 
and WZ terms vanish. The f n = 93 MeV is the pion 
decay constant. Here the space-time-dependent matrix 
field U(r,t) G SU(3) takes the form 



U(r,t) = A(t)U(f)A^(t), 



(5) 



where U(f) is the SU(3) matrix in which the Skyrme 
SU(2) ansatz is embedded: 



U(f) 



exp(iT-nF(r)) 




(6) 



The time-dependent collective coordinate matrix A(t) € 
SU(3) defines the generalized velocities A'(t)A(t) = 
t; Y] a —i A Q d a and the profile function F(r) is interpreted 
as a chiral angle that parametrizes the soliton. 



In this work we use the arctan ansatz for F(r) |13|: 



F(r) = 2 arctan (ro/r) 



(7) 



Here ro - the soliton size - is the variational parameter 
and the second power of ro/r is determined by the long- 
distance behavior of the massless equations of motion. 
After rescaling x = ref^, one obtains ro/r = xq/x. The 
quantity xq has the meaning of a dimensionless size of 
a soliton and it is determined by minimizing the classi- 
cal mass E c i. All relevant integrals involving the profile 
function turn into an integral representation of the Euler 
beta functions, which can be evaluated analytically . The 
accuracy of this method with respect to the numerical 
calculations is of the order of a few percent. In the chiral 
limit of the SU(2) Skyrme model, we obtain xq — \/T5/4 
and the arctan ansatz reproduces nucleon static prop- 
erties well [00] • Moreover we gain an insight on how 
different quantities depend on the soliton size, which in 
turn is a function of the symmetry breaker and e. 

In the SU(3) extended Lagrangian (||) we have a new 
set of parameters, namely, x = 36.4, (3 1 = —2.98 x 
l(T 5 GeV 2 , 8' = 4.16 x lCT 5 GeV 4 , determined from the 
masses and decay constants of the pseudoscalar mesons 
H . Owing to the presence of the /3' and 5' terms in E c i , 
xq becomes a function of e, / w , /3', and 5' , and it is equal 



„'2 
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where we use the symbol x' to distinguish it from the 
SU(2) case. After introducing the SB terms into the La- 
jangian (|J) , one can either treat them as a perturbation 
or one can try to sum up the perturbation series by 
numerically diagonalizing the resulting Hamiltonian 0]. 

The fitting procedure employed in this work is based 
on taking the physical values for f v and fx which takes 
care of the SU(3) symmetry breaking. In fact the SB af- 
fects the calculations in three different ways: (a) through 
the soliton size x' ; (b) via the explicit SB in the currents; 
(c) through the admixture of the higher SU(3) represen- 
tations in the baryon wave functions. It was shown in 
Rcf. [[ITJ] that the latter contributions to NHD are small, 
since the higher representations enter with small weights. 
Our estimate shows that they are of the order of 15 %. 
This uncertainty is of the order of the accuracy of the 
model which is reflected in the variation of the Skyrme 
parameter e depending on which static property is used 
in the fitting procedure. In the remainder of this paper 
we use the SU(3) symmetric wave functions. 
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For the evaluation of NHD, the important baryon 
static properties are the octet-decuplet mass splitting A 
and the axial decay coupling constant gA- The value of 
the only free parameter e sa 4 was successfully adjusted 
to the mass difference A of the low- lying l/2 + and 3/2 + 
baryons j^j. However, if we fix A, the constant gA is 
underestimated. This is a well-known problem of the 
Skyrme model, which can be cured in the more sophisti- 
cated chiral models involving quarks pl| . 

Therefore, we determine two values of the charge e 
through fixing A and gA to their experimental values. 
The arctan ansatz gives 



A = 



2A c (;Eq) 

14tt / /2 \ 
9A = 773 [2x o + 



15e 2 



n 16tt/3' , 2 7V2N, 
(1 — x) „ .„ x 



2 f 2 X 



where 



A 8 (aco) = 



225e 2 / 2 



i>(l + 2^| .,- 



4 ( 1-2(1 + 24)^ | 



192e/»A s (iJ)' 



25 , 



/3 



(9) 
(10) 

(11) 
(12) 



The quantity A c (x ) represents the rotation moment of 
inertia in coordinate space, while the A s (xg) is the mo- 
ment of inertia for flavor rotations in the direction of the 
strange degrees of freedom, except for the eighth direc- 
tion Hf7j . The static kaon fluctuations were omitted |h| 
in the derivations of equations (|| - |l^) . 

For the Lagrangian C, we calculate the matrix element 
of the product of two (V — A) currents between the octet 
states using of the Clebsch- Gordon decomposition M: 



< B 2 \d {SK) \B 1 >=<S> SK xJ^Cr 



(13) 



where <& SK is a dynamical constant and Cr denote the 
pertinent sum of the SU(3) Clebsch-Gordon coefficients 
in the intermediate representation R. The total matrix 



element is simply a sum < 6\ SK ' + 0\ * ' + Q\"~' >, 
with i = 1,...,4. The quantities $ are given by the 
overlap integrals of the profile function. Using the arc- 
tan anzatz (|7]), we obtain analytical expressions for the 
integrals as functions of x' : 



UWZ) 



■.(SB) 



~ 512 7T 2 x' 3[ U> ' 



4> 



SB 



Att 



V2 



64 xfi 



45 \ U 



(14) 



For the 0\ operator, R = 8 a , s or 27. ; then 



<Pt|6 1 |E +T >= 



wz 



+ 675 1 
1 / 7 , 



SK 



4 V 75' 



4> 



17 
1050 



(15) 



27 



4> 



SB 



The 27-piece is very small, which is an important proof 
of the octet dominance. 

By fixing A and gA to their experimental values, we 
obtain e = 4.228 and e = 3.385, respectively. In further 
calculations of the NHD amplitudes, we use the mean 
value e = 3.81 @ and x | e=3 . 8 i = 0.8782, i.e., 10% less 
than in the massless case. For f n = 93 MeV and e = 3.81, 
we obtain the following numerical values of the integrals 
@ in units of GeV 3 : 

<S> SK = 0.264, <f> wz = 0.004, $ ss = 0.005. (16) 



From eqs. @ and @ we find the following structure 
for a typical matrix element: 

<pT|Oi|s+t>= 

(-20.37$ SiC - 16.67$ wz - 27.38$ 5B )10~ 3 = 
(-5.38| S jc - 0.06|,yz - 0.14| SB )10- 3 Ge1/ 3 . (17) 



It is clear that on top of the octet dominance we also find 
the dominance of the Skyrme Lagrangian currents over 
the WZ and SB currents in the evaluation of a typical 
weak matrix element between two hyperon states. For 
e « 4, the SB and WZ terms are of comparable size and 
their coherent contribution to ( fL7| ) is below 4%. We see 
therefore from (|14|) that within this accuracy the result 



for (17) scales like 1/e (up to 5% due to the small vari- 
ations of the soliton size which weakly depends on e). 
The change of e between 3.385 and 4.228 produces 14% 
variations of the amplitudes ^4 C (0) and B v (rri^) around 
their mean values given in Table § 

In this work we have added factorizable, A s (rn^) and 
B (m 2 ), contributions to the Skyrme model amplitudes 
^4 C (0) and B p (mJ). The complete results are given in 
Table 0. Comparison of the total amplitudes A{rri^) and 
S(m 2 ) with experiment shows the following: 

(a) Short-distance corrections to the effective weak 
Hamiltonian are beyond doubt very important. 

(b) Signs and order of magnitudes of all amplitudes are 
always correctly reproduced. 

(c) s waves are in good agreement with experiment. 

(d) The Pati-Woo theorem violation Q and the 27- 
contaminations are found to be small. It is clear that 
the nonvanishing .A(E^) amplitude is still too small, in 
good accord with small values of the 27-contamination 
p3f , and that additional contributions are needed Jl2| . 

(e) p waves are subject to some uncertainties. Namely, 
in the Ref. |l4j] was shown that, in the Skyrme model, 
a contact term appeared and should be added to the re- 
sults for p-waves. That has been taken care of in Ref. 
p5[ and is not present in our approach. In our opinion, 
Z?(to 2 ) amplitudes are not fully described by our formu- 
las; nevertheless, they agree with experiment reasonably 
well. 
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(f) Finally, the factorizable contributions are small, 
and represent the fine tuning to the total amplitudes. 

To conclude, we would like to emphasize the fact that 
the pure Skyrme model Lagrangian C cannot explain 
nonleptonic hyperon decays pq| . However, the QCD- 
corrected weak Hamiltonian H^j' , together with the in- 
clusion of other possible types of contribution to the total 
amplitudes (K, X*-poles, and/or factorization; (l/2 ± *)- 
poles, etc.) supplemented by the Skyrme model, leads 
to a correct answer. This includes the explanation of 
the octet dominance, the |A/| = 1/2 selection rule, 
■A(£j) + 0, and the p/ s-wave puzzle. Nevertheless, this 
is certainly a matter of another series of studies. 
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TABLE I. The s-wave (A) and p-wave (B) NHD ampli- 
tudes. Choices (off, on) correspond to the amplitudes without 
and with inclusion of short-distance corrections, respectively. 
For the sake of comparison, we have added the constituent 
quark-model evaluation of the A c and B v amplitudes [9,10]. 



Amplitude(10~ 7 ) 




(A°_) 


(HI) 


(S + ) 




A c (0) 




off 


2.02 


-2.94 


-2.28 


0.02 






on 


3.84 


-5.56 


-4.34 


0.04 


A s {ml) 




off 


0.03 


-0.57 


-0.49 





§ 




on 


-0.42 


0.25 


-0.01 





A{ml) 




off 


2.05 


-3.51 


-2.77 


0.02 


(this work) 


on 


3.42 


-5.31 


-4.35 


0.04 


Exp. |2 






3.35 


-4.85 


-3.27 


0.13 


A c (0) 




off 


0.78 


-1.86 


-1.36 





CQM | 


1 


on 


1.49 


-3.53 


-2.59 





B ( p 1/2+) (m 


I) 


off 


20.1 


21.8 


13.7 


14.8 






on 


38.1 


41.4 


25.9 


28.2 


B s (ml) 




off 


3.6 


-1.5 


-0.4 





§ 




on 


6.0 


-2.4 


0.4 





B(ml) 




off 


23.7 


20.3 


13.3 


14.8 


(this work) 


on 


43.4 


38.2 


26.3 


28.2 


Exp. [[22 






22.3 


17.4 


26.6 


42.2 


B (l/2+)( m 


I) 


off 


2.9 


7.8 


7.3 


10.4 


CQM | 


\ 


on 


5.6 


14.8 


13.9 


19.7 
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